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Abstract 

We argue that a pseudo-one-dimensional electron gas is magnetized when a voltage bias is applied with the Fermi level tuned 
to be in the energy gap generated by a spin-orbit interaction. The magnetization is an indication of spin-carrying currents due 
to the spin-orbit interaction. The origin of the magnetization, however, is essentially different from the "spin accumulation" in 
two-dimensional systems with spin-orbit interactions. 
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1. Introduction 

In the present Letter, we consider a quasi-one-dimensional 
electron gas with the Rashba or Dresselhaus spin-orbit 
interaction [1,2,3,4]. We argue that the system is magne- 
tized when a voltage bias is applied longitudinally (Fig. 1) 
with the chemical potential tuned to be in the energy gap 
generated by the spin-orbit interaction. The magnetiza- 
tion indicates spin-carrying currents due to the spin-orbit 
interaction. The origin of the magnetization, however, is 
entirely novel; it is essentially different from the origin 
of the "spin accumulation" [5,6,7,8,9] in two-dimensional 
systems with spin-orbit interactions. 

A spin-orbit interaction affects the dispersion of a con- 
ducting channel of a pseudo-one-dimensional system in two 
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Fig. 1. (color online) A quantum wire under a voltage gradient. 



Fig. 2. (color online) A schematic view of the dispersion relations of 
the upper and lower subbands of a channel. The right-going current 
has a chemical potential fi^ = — |e| Vl of the left contact, while the 
left-going current has a chemical potential /xr = — |e|Va of the right 
contact. 

ways (Fig. 2): first, the dispersion of up-spin electrons and 
that of down-spin electrons, respectively shift sideways in 
the opposite directions; next, the crossing of the disper- 
sions at k% — — opens up an energy gap (see, however, 
the discussions in §4). Under this dispersion, we consider a 
simple theoretical state called the non-equilibrium steady 
state [10,11,12]; the right-going current has the Fermi dis- 
tribution of the left contact while the left-going current has 
the Fermi distribution of the right contact and they run bal- 
listically and independently without reflection. When the 
chemical potentials (the Fermi levels) of the right and left 
contacts are tuned to be in the energy gap of the disper- 
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sion (the case (a) in Fig. 2), the number of up-spin elec- 
trons in the left-going current and the number of down-spin 
electrons in the right-going current can differ, and thereby 
appears the magnetization. This argument does not apply 
when the chemical potential is in the middle of a subband 
(the case (b) in Fig. 2), where the right-going up-spin elec- 
trons cancel the magnetization of the right-going down-spin 
electrons. 

Magnetic properties of mcsoscopic systems with spin- 
orbit interactions are of great interest recently, particularly 
from the viewpoint of spintronics [13,14,15]; we could con- 
trol the dynamics of spins with an external electric field. 
Most studies so far have focused on two-dimensional sys- 
tems. The spin accumulation [5,6,7,8,9] and the spin-Hall 
effect [14,16,17,18,19] are two of the phenomena much 
discussed in two dimensions. The recent study by Eto et 
al. [15] showed the occurrence of a spin-polarized cur- 
rent in a pseudo-one-dimcnsional system with the Rashba 
spin-orbit interaction. This let us realize that pseudo-one- 
dimensional systems also contain much interesting physics 
of spin-orbit interactions. 

We may say that the origin of exotic phenomena due 
to spin-orbit interactions in pseudo-one-dimensional sys- 
tems is essentially different from the one in two-dimensional 
systems. Both the spin-polarized current predicted by Eto 
et al. [15] and the magnetization predicted here are due 
to an avoided crossing in the dispersion of the pseudo- 
one- dimensional system. We also note that we can treat 
pseudo-one-dimensional systems much simpler than two- 
dimensional systems. In fact, the linear response of the 
spin Hall effect [14,16] is recently under debate; it was 
shown [17,18] that the spin Hall coefficient vanishes in the 
DC limit ijj — > of the electric field under the presence of 
impurity scattering. It was then argued recently [19] that 
the long-range impurities may not eliminate the spin Hall 
coefficient even in the DC limit. In contrast, the treatment 
of the electronic transport in one dimension is well estab- 
lished. Our argument here is based on a common method 
widely used in studies of electronic transport in one di- 
mension [20]. We have already confirmed that inclusion 
of elastic impurities does not alter our conclusion qualita- 
tively [21]. 

2. Analytic calculation for a continuum model 

Let us begin our explanation with the Rashba system. 
The treatment of the Dresselhaus system is not much dif- 
ferent. The Hamiltonian of the system with the Rashba 
spin-orbit interaction is given by [22] 

p 2 + p 2 

U = 2m* " + aRSO ^ x ° y ~ ^ y&x ^ ' ^ 

where m* denotes the effective mass of an electron and 
c*RSO denotes the strength of the Rashba interaction. From 
a simplistic point of view, we could say that the term p y a x 
generates the gap at k x = while the term p x a y helps the 



current in the x direction to induce the magnetization in 
the y direction. 

We choose the y direction as the quantization axis of the 
electron spin, and hence use the following representations 
hereafter: 

*-(")• * ,= (L°i)' ,2) 

We assume that the system has only one mode in the z 
direction. The system is of length L in the x direction and of 
width W in the y direction with W <C L <C mean free path. 
Hence the electrons run along the wire ballistically through 
a few channels. 

We first diagonalize the Hamiltonian (1) in the momen- 
tum space with the bases 

\k x ,k y ,v v ) = -j±=e« k -* + W\v v ), (3) 

where k y — 0, ±2tt/W, ±4tt/W, . . . and a y = j, j. (For sim- 
plicity, we here imposed periodic boundary conditions in 
the y direction, which might be indeed materialized as a 
nanotube [23]. See the discussions in §4 below for fixed 
boundary conditions.) The Hamiltonian (1) is given by a 
two-by-two matrix in the spin space in the form 

ft 2 

H = - — - [k 2 + 9k (a y cos <j> - a x sin <p)] , (4) 

Zi 771 

where k 2 = k 2 + fc y 2 , (k x ,k y ) — fc(cos sin <p) and = 
2m* ttRso / fr- The spin rotation U = exp(z0<r z /2) gives the 
representation 

UnU^ = ^(k 2 + 0ka v ) (5) 
with the eigenvalues [24,25] 

e±{k x ,k y ) = ^{k 2 ±6k), (6) 

where the positive sign denotes the upper subband and the 
negative sign denotes the lower subband. For the lowest 
channel k y — 0, the dispersion has two branches of parabo- 
las with up and down spins, respectively. In the higher 
channels k y > 0, the two branches are mixed around the 
avoided crossing at k x = (Fig. 3(a)). Note that, away 
from the avoided crossing, the lower branch in the region 
k x 3> still predominantly has down spins while it almost 
has up spins in the region k x <C 0. 

Incidentally, lower channels arc nearly degenerate in the 
present system as shown in Fig. 3(a). This near degeneracy 
may be lifted by a confining potential U(y). As another 
remark, the double minima of the lower branch shown in 
Fig. 3(a) vanish in higher channels and the dispersion shifts 
upwards in even higher channels (Fig. 3(b)). The algebra 
hereafter slightly changes in the latter case but the final 
result (12) below is still valid. 

The right-going current originated in the left contact 
with the chemical potential contains all the states with 
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Fig. 3. (color online) Schematic views of the dispersion relation (6) 
for (a) some of the lowest channels and for (b) higher channels. The 
lowest channel has a real crossing at k x = 0, but a higher channel has 
an avoided crossing. In (b), the double minima of the lower branch 
vanish for a high channel and the whole dispersion shifts upwards 
as we go to even higher channels. 

positive group velocities. The group velocity in the x direc- 
tion is 



_ 1 de±(k x ,k y ) 
h dk x 



(7) 



which is positive wherever the slope of the dispersion is 
positive (the solid curves in Fig. 2). Each state (k x , k y ,±) 
has the magnetization per unit area in the y direction of 
the form 



m±(k Xl k y ) = ± 



Me 

LW 



■ cos 6 = ± 



Mc k x 



(8) 



because the spin rotation U changes the spin operator a y 
to o y cos 4> + &x sin <j). Thus we have the magnetization per 
unit area of the right-going current as 



m n (k y ) 



+ / m_/( £ _;T, M L) 



\-fc_ min fc_ r 



2tt/L 



oo 

j m+ f{e + ;T : ^)^- L 




(9) 



for a channel, where f(e;T,fi) denotes the Fermi distribu- 
tion with the temperature T and the chemical potential /i 
and we defined for each channel the following variables: 



k — r 



X {ky) = ^/A 



2/)2 



8m* 



£±o{k y ) 



2m* 



(k v 2 ±ek v ) 



(10) 



see Fig. 2. In changing the integration variable from k x to 
e in Eq. (9), we used the density of states D± given by 



1 



2tt 

T 



8ea 



dk x 



nfr \k x 
m*L k 



2 



r, 2 



-£±. 



(11) 



The van Hove singularity is canceled in Eq. (9) by the zero 
of m±. We can likewise obtain the magnetization m^ky) 
of the left-going current. The total magnetization per unit 
area is given by 



M t ot = ^ m R (fc a ) + m L (k y 



e+o(k y 



Me V- 

InhW 



/(e;l> L )-/( £ ;T,/z R ) 



y/uRso 2 + 2e/m* 



ds. 



(12) 



Note that only the integration over the energy gap survives. 
This is consistent with what we described in Fig. 2. 

Let us analyze the linear response. We define the 
chemical-potential bias as 



Ml = M + — and MR = M j~ 



(13) 



with AV = Vl- Vr = -A/x/|e|. The expansion of Eq. (12) 
with respect to A/j, is followed by 



X + (ky) 



ky ky 

where 

B(ky 

with 



(14) 



x-o(ky) 



1 



\/anso 2 + 2(k B Tx + fi)/m 

y '~ Av, T 



dx (15) 



and x ±0 (k y ) = £± ^ V i M -(16) 



[2 cosh(x/2)] 2 
3. Numerical demonstration 



Figure 4 shows the result of numerical calculation, where 
we used the values for an InGaAs/InAlAs heterojunc- 
tion [26]: a R sofr = 3 x 10" n [eV m] and m* = 0.041m e . 
We also set W — l[/im]. We find in Fig. 4(a) a peak of the 
magnetization around the energy gaps of the dispersion 
relation. We observe for AV = 1 [mV] the magnetization of 
the order of lfr 10 [J/T/m 2 ], which translates to 1(T 5 [G] 
if we assume that the thickness of the hetero junction is 
10[nm] [26]. 

Figure 4(b) shows the contribution of each channel, 
M„ = M ( - 1 \2mr/W), to the total magnetization M to t at 
T = 1 [K] . We note that each contribution is finite over the 
range of the energy gap of the respective channel, which 
is indicated by each bar. In the present case, some of the 
energy gaps overlap and hence a sharp peak in Fig. 4(a). 
Figure 4(b) also tells us that the ripples in the tail of the 
curve for T = 1[K] in Fig. 4(a) is due to the summation 
over the channels. 

Going back to Eq. (12), we show in Fig. 5 the dependence 
of the magnetization on the voltage bias AI^ = — |e|A/i be- 
yond the linear response. The magnetization grows mono- 
tonically in this regime. Although the contribution from a 
channel saturates as the voltage bias exceeds the energy gap 
of the channel, more and more channels contribute to the 
summation in Eq. (12) and hence the monotonic increase. 
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Fig. 5. (color online) The magnetization dependence on the voltage 
bias AV for fi = 0[meV]. 



10 

5 

I » 

-5 



-10 



-3 -2.95 -2.9 -2.85 -2.! 



Fig. 4. (color online) (a) The response of the magnetization per unit 
area to the voltage bias AV. (b) Contributions M n = M' 1 ' (2nn/W) 
from various channels for T = 1[K]. The bars indicate the energy 
gaps of the channels. In higher channels, the double minima at 
k x = fc_ m i n vanish and the "energy gap" shifts to the right. 



4. Numerical calculation for a ladder model 

So far, we have assumed periodic boundary conditions in 
the y direction. When we introduce a confining potential 
U(y) into the Hamiltonian (1), the gap at k x = may be 
closed; the standing wave in the y direction includes both 
positive and negative k y and thus the gap-opening effect of 
the term p y a x is cancelled out. We can then revive the gap 
by introducing a magnetic field in the x direction: —H x a x . 
In order to confirm our prediction in the confined geometry, 
we carried out a numerical calculation of the corresponding 
tight- binding model on a two-leg ladder [27]: 



«=-* E E E 

x=-oo tr=t,4, Vy=l,2 



Fig. 6. (color online) The magnetization in the y direction of the 
model (17) for t = 1, da = 0.2 and H x = 0.01. There are two 
subbands in this region, just as in Fig. 2. The broken curves in the 
negative side indicates the contribution of the right-going states with 
k x > in the lower subband (the second integral of Eq. (9)), the 
broken curve in the positive side indicates the contribution of the 
left-going states with k x > in the lower subband (the first integral 
of Eq. (9)), and the dotted curve indicates the contribution of the 
left-going states in the upper subband (the third integral of Eq. (9)). 
The solid curve is the sum of the three contributions, which survives 
only in the energy gap. 

+ c x . 2 Je- l9a ^c xA , a + K.c^j 

oo 

x= — oo y=l,2 <T=t,J, 

where a is the lattice constant. Figure 6 shows a result for 
T = 0, confirming our prediction. Note here that the width 
of the energy gap is of the order of H x but the magnitude 
of the magnetization is of the order of 8. Details of the cal- 
culations on the lattice system will be reported elsewhere. 

The magnetic field H x applied in the present section in 
order to revive the energy gap at k x = may induce ad- 
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ditional spin polarization at the left and right ends of the 
system, namely at the contacts to the baths. Even without 
the magnetic held, the dissipation at the contacts might 
generate some form of spin polarization. The precise mag- 
nitude of the polarization would strongly depend on the in- 
teraction between the system and the baths and is beyond 
the scope of the present paper. The magnitude, however, 
should be of the order of 0(i _1 ) compared to the magneti- 
zation Mtotj the magnetization given in Eq. (12) and Fig. 6 
is the one per unit area, that is, a bulk quantity. Wc hence 
assume here that the spin polarization localized around the 
contacts is negligible. 

Finally, the Hamiltonian with the Dresselhaus interac- 
tion takes the form [22] 



H = 



Px + Py 

2m* 



«DSO (Px&x - PyOy) 



(18) 



The rest of the formulation is the same as above with 9 = 
2m*aDSo/fo and with the spin quantization in the x direc- 
tion as in 



1 

-1 



-i 

1 




The magnetization in the x direction is measured but the 
expression is the same as Eq. (12). 

5. Summary 

To summarize, wc predicted that the magnetization ap- 
pears under a voltage bias in pseudo-one-dimensional sys- 
tems when the Fermi levels are tuned to be in the energy 
gap due to spin-orbit interactions. The magnetization is an 
indication of the spin-carrying current, but its origin is es- 
sentially different from the one in two-dimensional systems 
with spin-orbit interactions. The energy gap due to spin- 
orbit interactions plays a key role in exotic phenomena of 
pseudo-one-dimensional systems. Without taking account 
of scattering processes, the argument may be naive; we nev- 
ertheless believe that it is worth reporting and should be 
checked experimentally. We can also argue that the mag- 
netization appears under a temperature gradient [28,29,30] 
rather than a voltage gradient; this effect will be reported 
elsewhere. 
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